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We investigate the effects of the dynamic bosonic fluctuations on the Fermi surface reconstruction 
in two dimensions as a model for the underdoped cuprates. At energies larger than the boson energy 
uJb, the dynamic nature of the fluctuations is not important and the quasi-particle dispersion exhibits 
the shadow feature like that induced by a static long range order. At lower energies, however, the 
shadow feature is pushed away by the finite u)t>- The detailed low energy features are determined by 
the bare dispersion and the coupling of quasi-particles to the dynamic fluctuations. We present how 
these factors reconstruct the Fermi surface to produce the Fermi arcs or the Fermi pockets, or their 
coexistence. Our principal result is that the dynamic nature of the fluctuations, without invoking a 
yet-to-be-established translational symmetry breaking hidden order, can produce the Fermi pocket 
centered away from the (n/2,n/2) towards the zone center which may coexist with the Fermi arcs. 
This is discussed in comparison with the experimental observations. 
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I. INTRODUCTION 

The "Fermi arc" picture was advanced by the angle- 
resolved photo-emission spectroscopy (ARPES) to under- 
stand the enigmatic pseudo-gap state in the underdoped 
cuprates.— ~— The ARPES, with its momentum resolution 
capability, established that in this pseudo-gap state the 
gapped region is mainly in the (0, 7r) and (%, 0) region 
while the Fermi surface (FS) exists in the diagonal direc- 
tion. Then, the picturesque view of the pseudo-gap state 
is that the gapless portion of the FS forms an open ended 
arc, rather than a closed loop as in ordinary metals. It is 
extremely difficult to understand the abrupt truncation 
of the FS in the Brillouin zone. The Fermi arc has thus 
puzzled the physics community and triggered enormous 
research efforts £ 

This Fermi arc picture was challenged by the observa- 
tions of the quantum oscillation under the applied mag- 
netic field H The transport and thermodynamic prop- 
erties exhibit the periodic oscillations as a function of the 
inverse magnetic field. The standard interpretation is in 
terms of the closed loop of the FS, or, the Fermi pockets. 
The oscillation is due to the quantizated Landau levels, 
and its periodicity is proportional to the area of the Fermi 
pocket. It is found to be only a few percent of the FS 
area of optimally or overdoped cuprates. In the theory of 
usual metals, such a small FS would require a change of 
translational symmetry from overdoped to underdoped 
cuprates. The problem is that there is no direct experi- 
mental evidence for the translational symmetry breaking 
for the compounds exhibiting the small FS. Moreover, the 
Fermi pocket is at odds with the Fermi arc picture from 
ARPES. Although the ARPES were done above T c with 



no magnetic field and the quantum oscillations in the 
low T and strong external field, the views they advance, 
the Fermi arc and Fermi pocket, seem contradictory each 
other and need to be reexamined. 

The recent laser ARPES on the single layer 
P^Si^-jTa^CuOg compounds by Meng et al£ is indeed 
very interesting in this regard. They observed, with the 
improved resolution, that the ungapped portion of FS 
forms a closed loop, e.g. the Fermi pocket, rather than 
the Fermi arcs at the doping levels of 11 and 12 % for 
Bi2Sr2- 2; La a ;Cu06 . Moreover, the center of the Fermi 
pocket is shifted from the (tt/2, 7t/2) toward the zone cen- 
ter (r point). The translational symmetry breaking, let 
alone its yet-to-be-established existence, can not explain 
their results because a salient feature of the reconstructed 
FS induced by the broken translational symmetry of pe- 
riod doubling is that the FS is symmetric with respect to 
the (tt, 0) - (0, tt) line. 

Here, we wish to understand the Fermi pocket centered 
away from the (7r/2,7r/2) point without invoking the 
translation symmetry breaking in terms of the dynamic 
bosonic fluctuations. We first consider a dynamical col- 
lective mode coupled with quasiparticles at the antiferro- 
magnetic (AF) wave vectors only (the correlation length 
£ — > oo) for simplicity and illustration of basic ideas. 
Then, the more realistic cases of finite £ are presented 
with self-consistent numerical calculations. 

There have been many attempts to understand the 
Fermi arcs and pockets in the cuprates. Each of them has 
discrepancies with the experimental observations such as 
the shape, location, or the spectral weight*^— On the 
other hand, the dynamic nature of the bosonic fluctua- 
toins peaked at (7r, 7r), without invoking a hidden order 
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which breaks the translational symmetry, can produce 
the FS evolution from the large FS to Fermi arc to Fermi 
pocket as the coupling is increased. More specifically, 
it can induce (1) the Fermi pocket centered away from 
the (7r/2, 7r/2) towards the T point, (2) the ratio of the 
spectral weight at the back side of the Fermi pocket to 
the inner side is about 10 -2 , (3) coexistence of the Fermi 
pocket and the large main FS, and (4) the dispersion kink 
along the nodal direction at energy w 0.05 eV. These 
are in agreement with the recent laser ARPES experi- 
ment of Meng et al- and numerous previous experimental 
reports*^— 

After the bare band dispersion is determined there are 
three factors which affect the Fermi surface reconstruc- 
tions: the fluctuations correlation length £, coupling con- 
stant a, and the boson frequency scale ujb- More discus- 
sion about their possible microscopic origin and relation 
will made later in Sec. V in connection with other ap- 
proaches. For now, we first take the Einstein mode of 
Wf, for simplicity, Ub = 0.05 eV was chosen to match 
the kink energy»i2r— We will also consider the realis- 
tic frequency dependent bosonic spectrum recently de- 
duced by Bok et a& by inverting the laser ARPES on 
Bi 2 Sr2CaCu208+,5 . We then perform detailed numeri- 
cal calculations and show that the dynamic nature of the 
collective mode can account for the FS evolution without 
introducing a yet-to-be-established hidden order param- 
eter. 



II. IDEA AND FORMULATION 

We consider the renormalization of the fermions due to 
the coupling to the dynamic bosonic fluctuations F(q, w) 
with the coupling vertex a(k, k'). The self-energy of the 
fermion is given by2i 
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where A is the spectral function of the fermion, and / 
and n are the Fermi and Bose distribution functions, re- 
spectively. 



A(k',e) = --Im- ' 
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We took the fluctuation spectrum of the following factor- 
ized forms^i 



E 



Q— ±7r/a,±7r/a 



a 2 F(k,k',e') = a(k,k') 2 i ;l (e' 
T/tt T/tt 

(fe-Q,) 2 + r 2 fe-Q^ + r 2 



, (4) 



where a is the lattice constant, q = k' k, and V = 7r/£. 
The coupling a may depend on the wave vectors k and 
k', but for simplicity we will consider a constant a and 
the Einstein model of frequency ujb first. 



a 2 F(e') = a 2 [S(e' - u b ) - 5(e' + u b )] 



(5) 



Some remarks will be made on the more realistic fre- 
quency dependence of -F(q, e') and the momentum de- 
pendence of a(k,k') later. Eqs. ([T]) and ([2]) constitute 
the coupled self-consistency equations. They are solved 
self-consistently for the self-energy via numerical itera- 
tions. A very similar problem was investigated by Grilli 
et al. for the one dimensional electronic systems *22, It is 
extended to two dimensions in the present work fully self- 
consistently. 

Let us first consider the simple case of T — > and T — > 
to gain underlying physics. That is, the boson mode 
is of a delta function in both the energy and momentum 
channels. Then, in the limit T — > 0, Eq. (p} is reduced 
to 
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where kg = k + Q and & is the step function. A useful 
approximation is to take 



We then have 

E(k,w 
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ik Q = tk Q + sgn(£ kQ )uj b . (10) 
The Green's function of quasi-particle (qp) is given by 

1 
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This form of the Green's function appeared previously in 
the context of the pseudogapi 12 i 14 The coupling vertex 
a(k, k') of present approach corresponds to the pseudo- 
gap Ar of Ref.— . It will be interesting to check to what 
extent this mapping is valid. An important distinction of 
the present approach is that the dynamics of the bosonic 
fluctuations is explicitly built in via ^ Q of Eq. (TTU|) . It 
is precisely this dynamics which gives rise to the Fermi 
arcs as we will see now. 
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The qp dispersion E(k) is determined by 

G- 1 (k,uj)=uj-^ k -^(k,uj) = 0, (12) 
which gives 

1 



E ± (k) = 



(13) 



The results may approximately be extended to the case 
of finite correlation length l/£ ^ following Ref.— by 
replacing the imaginary part of the frequency by <5 = 

The Green's function may be cast into the form 



G(k,cj) 



uj + iS — E + uj + iS — E- 
where the coherence factors are given by 



(14) 
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The -E+ and E— represent, respectively, the electron and 
hole bands. The spectral function A(k, uj) is then 

A(k,u>) = —ImG(k,u) 

= ulS(u> - E + (kj) + v 2 k S(uj - E_{k)). (16) 

The spectral function is directly probed by the ARPES. 

III. PRELIMINARY ANALYSIS 

Before showing the detailed numerical results, we will 
first present the preliminary analysis to gain underly- 
ing physics of the problem. The bare dispersion of 
F^S^-zLa^CuOg is taken as 



£fc = — 2t[cos(k x a) + cos(k y a)} + At' cos(k x a) cos(k y a) 
— 2t" [cos(2k x a) + cos(2fc a a)] — fj,, 



(17) 



where t = 0.25, = 0.058, t" = t'/2 eVU The FS 
corresponding to the £/. and E±(k) with a = 0.1 eV, 
u!b = 0.05 eV, and /i = —0.208 eV corresponding to the 
slight underdoping of 12 % are shown in Fig. 1. The 
nodal cut of k x = k y and several cuts parallel to it are 
also shown with dashed lines. 

Along the cuts the qp dispersions are presented in Fig. 
2 to better reveal the dynamically generated gap close 
to the shadow FS. Fig. 2(a) is the hole band dispersion 
_E_(k) in solid blue and electron band E+(k) in dashed 
green lines along the nodal cut given by Eq. (IT3"|) . The 
important point is that the hole band dispersion exhibits 
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FIG. 1: The reconstructed FSs for uj b = 0.05, a = 0.1, and 
fi = —0.208 eV corresponding to x — 12 % doping. The thick 
black curve is the bare FS determined by = of Eq. (|17[) 
and the green curve is the shadow FS of £k+g = 0. The 
red curve around the (tt/2,-k/2) point shows the hole FS of 
E-(k) — 0. Notice that the outer portion is gapped and the 
hole FS forms the Fermi arc. The blue curves around (0, tt) 
and (tt, 0) are the electron FS of E+(k) = 0. As a increases it 
may disappear as shown in Fig. 3. The nodal cut of k y — k x 
and the parallel cuts of k y — k x + O.ln with n = 1 — 7 are 
indicated with dashed green lines. The band dispersions along 
these cuts are shown in Fig. 2. 



the abrupt jump at k x a/i: « 0.6, or the gap of about 
2u!b- The dynamically induced gap was noticed by Grilli 
et al. for the one dimensional electronic systems^ 

The gap of 2u>b means that for uj = there exists only 
a single k point which satisfies w = E-(k), while for 
\ui\ > tub there exist two k points, one close to the orig- 
inal FS and the other to the shadow FS. That is, the 
shadow feature is present for \u\ > w&, but is absent for 
uj = 0. This is in accord with general expectations: A 
physical system may have long (but finite) ranged order- 
parameter spatial correlations which fluctuate with the 
frequency u>b- The system then appears to be ordered 
above 0Jb- For energies larger than u>b with respect to 
the Fermi energy the spectra should resemble an ordered 
system. On the other hand, at lower energies electrons 
"sense" averaged order-parameter fluctuations, and the 
system appears to be not disturbed much from the one 
without the collective mode. 

The blowup of the hole band dispersions is shown in 
Fig. 2(b) along the cuts parallel to the nodal cut. Notice 
that the gap survives beyond k y a/-K = k x a/-n + 0.4. It 
simply means that the gapless portion of the FS forms an 
open ended arc as shown with the thick red solid curve 
in Fig. 1. We stress that the abrupt truncation of the 
FS, which seemed so puzzling, is naturally understood in 
terms of the dynamic boson mode. 

Fig. 3 is the 3D plot of the spectral function A(k, uj) 
as a function of k at uj = 0. Fig. 3(a) is for ujb — 0.05 eV, 
a = 0.1 eV, and x — 12 %. Because of the dynamically 
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FIG. 2: The hole and electron band dispersions, E-(k) and 
£/+(k), along the cuts parallel to the nodal cut indicated with 
the green dashed lines in Fig. 1 with the same parameters. 
Fig. 2(a) shows the hole band in solid blue and electron band 
dispersion in dashed green curve along the nodal cut. Notice 
that the hole band dispersion has a gap of 2cJt at k x a/n ~ 0.6. 
The detailed hole band dispersions are plotted in Fig. (b) 
along the cuts of k y a/n = k x a/n + O.ln with n = — 6 from 
right to left. The gap survives beyond k y a/n — k x a/n + 0.4. 
It means that the gapless portion of the FS forms an open 
ended arc as shown with the thick red solid curve in Fig. 1. 



generated gap close to the shadow FS discussed above, 
the spectral peak shows up only over a part of the FS 
instead of a closed loop as in Fig. 3(b). Also the spec- 
tral peaks from the electron band show up around (0,tt) 
and (it, 0) for a weak a. Now, the FS may evolve to 
the Fermi pocket as the coupling a is increased. As an 
illustration, Fig. 3(b) is the 3D plot of the spectral func- 
tion for a — 0.2 eV with all other parameters fixed. The 
Fermi pocket is clearly formed. The peaks from the elec- 
tron band are substantially reduced. Physics behind the 
Fermi arc/Fermi pocket induced by the dynamic fluctua- 
tions is quite simple: The self-energy correction given by 
Eq. (P) dynamically generates a gap close to the shadow 
FS of magnitude of about 2u>b, marked by "2w£," in Fig. 
2(a). As a increases, the gap between the electron and 
hole bands marked with "2cv" in Fig. 2(a) becomes larger 
and the hole dispersion £L (k) of Eq. (|T3| is pushed down. 
Consequently, the qp states above the 2uib marked with 
"A" in Fig. 2(b) touch the FS. Then FS forms over a 
closed loop, which is the Fermi pocket. 

Where -E-(k) = either Fermi arc or Fermi pocket 
shows up. If two k points satisfy £L(k) = along any 
cut between the two hot spots and parallel to the nodal 
cut, then the Fermi pocket is produced. If, on the other 
hand, either one or two k points satisfy £L(k) = 0, then 
a portion of a pocket is missing, which is just the Fermi 




FIG. 3: 3D plot of the spectral function A(k,L) = 0). Fig. (a) 
is for Wi,=0.05, a = 0.1, and /i = —0.208 eV corresponding to 
x = 12 %. The Fermi arc appears because of the dynamically 
generated gap of magnitude of 2u?b close to the shadow FS. For 
Fig. (b), the parameters are the same as Fig. (a) except a = 
0.2. The Fermi pocket appears now because the gap is pushed 
down below the Fermi energy. Fig. (c) is for |a(k, k')| 2 = 
QfQ |k x k'| with ao = 0.3 and all other parameters the same 
as (b). The formation of the Fermi pocket is less favored. 



arc. Both cases can be produced with the simple formula 
of Eq. ([Dp]) depending on the parameters as discussed 
above. 

Also interesting is the relative weights of the two peaks 
of the Fermi pocket. For example, along the nodal cut, 
there appear two peaks near the main band and shadow 
band as a function of the momentum amplitude. The 
ratio of the spectral weight on the back side of the pocket 



5 



to that on the main FS is from Eq. (|T5 



«2(& = o) 



A'i,> 



0.01 



(18) 



in accord with the experimental observation. 9 

We also considered the momentum dependent coupling 
a(k, k') as suggested by Varma and coworker o 25 i 26 and 
also by Yang et alH 



|a(k,k 



/\|2 



al |k x k'| 



(19) 



This form of coupling will modify the qp disperion less 
along the nodal cut because k x k' k there. The Fermi 
pocket formation is less favored. In Fig. 3(c), we show 
A(k, uj — 0) for ao = 0.3 eV with all other parameters 
the same as Fig. 3(b). The Fermi arcs are formed instead 
of the Fermi pocket as anticipated. 

From the shapes of the Fermi arcs shown in Figs. 1 
and 3, one may notice that the arcs turn in near the 
ends. It means that the Fermi arcs seem to deviate from 
the underlying FS near the ends. Norman et al. argued 
that this is a generic feature of the pseudogap induced by 
q ^ order parameters^ This point seems to apply to 
the Fermi arcs induced by dynamic fluctuations as well 
although the turning in looks weaker. 

Now we understood the basic physics underlying the 
Fermi arc and Fermi pocket formation with the simple 
dynamic bosonic fluctuations of Ub = 0.05 eV and £ = 
oo. But, as £ — > oc the boson mode must get soft and 
approach ljj, — > 0. This relation was not satisfied in the 
simple case just presented. We therefore performed the 
full self-consistent calculations in the following section 
with finite £ and temperature. The important message of 
the numerical calculations will be that the dynamically 
generated gap of 2w& in the back-side of the pocket as 
shown in Fig. 2(a) remains intact as can be seen from 
Fig. 4(a). It means that the qualitative feature of the 
FS evolution from the large FS to Fermi arcs to Fermi 
pockets is unaffected. This is easy to understand. The 
magnitude of discontinuity being determined by Ub , it is 
insensitive to £ = oo or not as presented in the following 
section. 



IV. NUMERICAL RESULTS 

The previous discussion is based on approximate solu- 
tion of the self-energy of Eq. (jHJ. Although the approx- 
imation permits the simple and useful results discussed 
in the previous section, some of the results may be an ar- 
tifact of the approximation. We therefore performed the 
full self-consistent calculations via numerical iterations of 
the coupled equations of Eqs. (p} and ©. We considered 
the finite correlation length (r = 7r/£ ^ in Eq. (QJ) and 
non-zero temperature. The more realistic frequency de- 
pendent a 2 F(e') as extracted by Bok et aL— is also con- 
sidered. The important effects of the self-consistency are 



that (a) the Fermi arc and Fermi pocket coexist and (b) 
the center of the Fermi pocket gets displaced towards the 
zone center. The fine details are determined by the pa- 
rameters like a, T, and T. The non-zero T, non-zero tem- 
perature, or the frequency distribution of a 2 F(e r ) smear 
the fine structures out. 

It is interesting to note that the laser ARPES exper- 
iments observe that the Fermi pocket coexist with the 
Fermi arc. The coexistence may be understood as fol- 
lows: Let us fist consider the hole FS. The electron FS 
follows the same arguments. The spectral function of Eq. 
(fT6f indicates that the peaks show up as a function of k 



where £L(k) = or v\ is maximum for the hole FS. The 
reconstructed hole FS may appear in the region where 
£fc > and £k+Q > around the (n/2,n/2) point. The 
loci of maximum v 2 can be seen most clearly in the limit 
a = 0. Inspection of the coherence factor v 2 of Eq. (fT5|) 
in the limit a = reveals that Vk — 1 for ^k+Q > £,k- 
Simultaneously, E- (k) needs to be close to as the delta 
function of Eq. (fTo) requires. Both conditions are satis- 
fied where = 0. It is expected that peaks are produced 
close to the original FS due to the v 2 factor of Eq. (f!5|) . 

The coexistence may also be understood as follows: 
The so-called two-pole approximation of Eq. ((5J produces 
two qp branches. Next order approximation, the three- 
pole approximation, is to use Eq. (|16[) to the self-energy. 
It produces three qp branches. Straightforward calcu- 
lations reveal that, along the nodal cut near u = for 
example, there exist one branch close to the bare FS, and 
two branches almost symmetric around the (jr/2, tt/2) at 
(tt/2 ± e,7r/2 ± e). Between the two, the one closer to 
the bare FS, (tt/2 — e,ir/2 — e), merges with the branch 
near the bare FS to form the main FS, and the one at 
(n/2 + e, n/2 + e) forms the back side of the Fermi pocket. 
The self-consistent calculations to be presented below 
maintain this feature to produce the coexisting Fermi 
arcs and pockets. 

Another effect of finite £ is to exhibit the dispersion 
kink near uj « — ujb- In the limit of £ — » 0, it is simple to 
see that 



S(k, uj) = a 2 In 



uib — w 



LUb + w 



(20) 



Then, the slope of the qp dispersion changes from 1 + 
2a 2 /ujb to 1 as w increases past u>b- This dispersion kink 
along the nodal cut was observed by many groups and 
has been the focus of intense debate. 

For finite £, the summation over k' in Eq. @ is not 
a delta function. The k' summation was performed by 
using the 2D FFT (fast Fourier transform) between the 
momentum and real spaces using the convolution relation 



£F(k'-k)G(k')=F(r)G(r) 



(21) 



2 8 points were taken for the FFT along each axis. For 
a not too large a convergence took about 10 iterations. 
For a larger than about 0.22 eV the procedure failed to 
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FIG. 4: The qp dispersion along the nodal cut from the self- 
consistent calculation, a = 0.18 eV, T — 0.02, and T — 
for Fig. (a) and T = 200 K for Fig. (b). In the ARPES 
experiments the features above uj = are cut off by the Fermi 
distribution function. Notice that the gap of 2w& opens up 
in the back side of the pocket for l/£ 7^ as well. Also 
notice the shadow feature around u) = induced by the finite 
temperature effects. The shadow band disperses away from 
the zone center in accord with the observation by Meng et al. 



converge in our numerical iterations. This could be an 
indication of a topological change of the Fermi surface. 

Fig. 4 is the density plot of the spectral function 
A(k, u>) along the nodal cut as a function of k x a/ir and 
w with a = 0.18 eV, T = 0.02, and T — for (a) and 
T = 200 K for (b). At T = the shadow band appears 
with the gap of 2wf, centered around the Fermi energy. 
The main band is modulated by the u>b and the gap of 2a 
is not distinguishable. The dispersion modulation, being 
determined by the energy Wf, in the case of the Einstein 
mode, is expected to be weakened if the spectrum has 
a finite energy distribution. This expectation is indeed 
the case as will be presented below in Fig. 7. Also note- 
worthy is that the shadow band disperses away from the 
(7r/2,7r/2) as the energy is lowered in accord with the 
ARPES observation of Meng et al. Compare with the 
lower row of the plots b-d of the Fig. 1 in the Ref A 

An important role of the finite temperature presented 
in Fig. 4(b) is to bring up the qp states below the Fermi 
energy (marked by "B" in Fig. 2(b) for the two pole 
approximation) to form the Fermi pocket. This is in 
contrast with the simple results presented in the previ- 
ous section. The non self-consistent preliminary analy- 
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FIG. 5: (a) The spectral function as a function of k x and k y at 
oj = with the same parameters as the Fig. 4(b). The Fermi 
pocket is formed because of the temperature induced shadow 
feature around w = 0. (b) The plots of the spectral function 
of (a) along the cuts parallel to the nodal cut. From right to 
left are the cuts of k y a/n = k x a/n + 0.2n with n = — 4. 



sis indicated that the qp states above the Fermi energy 
(marked by "A" in Fig. 2(b)) are pushed down by the 
a and form the Fermi pocket. This picture is modified 
in the self-consistent calculations: As a increases the qp 
dispersion above the Fermi energy bends back as can be 
seen from Fig. 4(a) to keep the gap as intact as possible 
because the total energy will be lowered by not occupying 
the higher lying states. Instead the shadow band disper- 
sion below the Fermi energy is extended above the Fermi 
energy to form a pocket as can be seen from Fig. 4(b). 
Note that at the Fermi energy the dispersion from be- 
low is closer to the zone center than the dispersion from 
above. Consequently, the pocket is displaced towards the 
zone center away from the (jv/2, it/ 2) point as shown in 
Fig. 5(a). 

In Fig. 5(a) we show the density plot of the spectral 
function A(k, u = 0) as a function of k with the same 
parameters as the Fig. 4(b). Note the formation of the 
pocket coexisting with the main Fermi surface. The cen- 
ter of the pocket is shifted to the zone center away from 
the (tt/2,it/2) point as discussed above. Fig. 5(b) is the 
plots of the spectral function of A(k, ui — 0) along the 
cuts parallel to the nodal cut. From right to left are the 
cuts of k y a/ir = k x a/ir + 0.2n with n = — 4. Note the 
small peaks near the back side of the pocket. The ratio 
of their spectral weights to those on the main bands is 
found to be about 10 -2 . 
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FIG. 6: The input Eliashberg function taken from Bok et al. 
The dimensionless coupling constant A m 1.5. It corresponds, 
for the Einstein mode of uit = 0.05 eV, to a ~ 0.19 eV. 

We now turn to the realistic frequency dependent fluc- 
tuation spectrum. It was taken from Bok et al£Q with 
a constant a. The input Eliashberg function a 2 F(uj) is 
shown in Fig. 6. The extracted fluctuation spectrum has 
a peak around u « 0.05 eV, flattens out above 0.1 eV 
and has a cut-off at approximately 0.35 eV. The dimen- 
sionless coupling constant 

f°° , 2a 2 F(uj) 

A = / duj (22) 

Jo w 

is A fa 1.5. The Eliashberg function a 2 F(9,u>) deduced 
by Bok et al., where 9 is the tilt angle with respect to 
the nodal cut and U) is the energy, is that the functions 
along different angles collapse onto a single curve below 
the angle dependent cut-off energy lu c (9). The cut-off is 
maximum along the nodal cut, uj c » 0.35-0.4 eV, and 
decreases as the angle is increased. In the present calcu- 
lations this angular dependence of the cut-off energy of 
the Eliashberg function was disregarded. 

Eqs. ([U, ([2]), and (|4]) were solved self-consistently via 
iterations taking the extracted a 2 F of Fig. 6 into con- 
sideration. The k' summation was performed using the 
2D FFT as explained above. The finite range of the fluc- 
tuation spectrum instead of a delta function is to smear 
out fine structures of the spectral function as can be seen 
by comparing Fig. 7(a) with the Fig. 4 of a delta func- 
tion fluctuation spectrum. In Fig. 7(a) we show the dis- 
persion along the nodal cut at T = 100 K, that is, the 
density plot of A(k,uj) as a function of k x a/ir and w. 
The shadow band is also smeared out and its width is 
increased as the energy is lowered. In Fig. (b) the den- 
sity plot of A(k, uj = 0) is shown as a function of k. The 
pocket becomes weaker compared with the delta function 
fluctuation spectrum of Fig. 5(a). 

In order is to make a comment on implication of the 
coupling constant A of Eq. (|22|) on superconductivity. 
The approximate T c formula for <i-wave superconductors 
is 

T c = u av e ~( 1+A ^, (23) 

where X s and Xd are the coupling constant in the normal 
and pairing channels, respectively. The A of Eq. (|2"21 is 



1 







FIG. 7: (a) The spectral function A(k, cj) as a function of 
k x a/ir and ui along the nodal cuts at T = 100 K. The Eliash- 
berg function of Fig. 6 was taken into consideration, (b) The 
3D plot of A(k,Ld) as a function of k at u = 0. 



A s because it was extracted in the normal state. T c ~ 
150 K is produced if we take A^ = 0.8A S . This is in 
accord with the expectation g = Xd/X s < 1 for d-wave 
superconductors^ 



V. REMARKS AND OUTLOOK 

We investigated the effects of the dynamic nature of 
bosonic fluctuations on the Fermi surface reconstruction 
as a model for the underdoped cuprates. The dynamic 
fluctuations induce the gap of magnitude 2w& close to the 
shadow Fermi surface as Fig. 2 demonstrates. Then, the 
Fermi surface in momentum space can be truncated un- 
like the Fermi surface reconstruction induced by a long 
range order. Therefore, the Fermi arcs are naturally 
induced by the dynamic fluctuations. The Fermi arc 
and/or Fermi pocket is formed as Figs. 3, 5, and 7 show 
depending on the coupling constant a or the tempera- 
ture T or the correlation length £. The Fermi pocket is 
formed by the filling in of the dynamically generated gap 
by the non-zero temperature or the energy distribution 
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of the bosonic spectrum a 2 F(w). The self-consistency 
enables the Fermi arcs and pockets coexist and moves 
their center towards the zone center. 

There have been many works along the same path 
adapted in this paper, namely, employing bosonic fluctu- 
ations to compute the renormalization of the electronic 
properties. See Refill for a recent review. Now, it will be 
in order to make some comments on and comparison with 
a few recent relevant works. In Ref,— , Greco computed 
the electronic polarizability of d density wave instability 
(or, flux phase) with the t-J model. It was used as the 
bosonic fluctuations to couple to the electrons. The cal- 
culation is non-self-consistent and assumes the true phase 
transition of the flux phase. The symmetry broken phase, 
however, is yet to be confirmed experimentally. Never- 
theless, Greco addressed some of the points we did not 
touch in this paper like the temperature dependence of 
the Fermi arc length.— In the absence of any symmetry 
broken phase in the pseudogap doping range, however, 
we did not specify the mechanism of the boson mode in 
this paper. Instead we took a phenomenological effective 
interaction between electrons like Fig. 6. Because our 
main point was to demonstrate the Fermi surface evo- 
lution with a, we did not pursue the questions like T*, 
temperature dependence of the arc length, and so on, 
leaving them as further studies. 

Dahm et al. made a check if a self-consistent descrip- 
tion is possible between ARPES and inelastic neutron 
scattering (INS) for YBa2Cu306.6'^ 1 In more detail, they 
fitted the INS to extract the spin susceptibility. Then 
they used it as the bosonic fluctuation to couple with the 
electrons to calculate the self-energy. The results were 
consistent with ARPES intensity and nodal dispersion, 
and the kink along the nodal cut was produced. This 
nodal kink is expected in their work because T of the 
extracted susceptibility is non-zero. 

The dynamic fluctuation model with no long-range or- 
der of the present paper successfully describes the FS 
evolution from the large FS to Fermi arc to Fermi pocket 
as the coupling is increased. Particularly, the enigmatic 
abrupt truncation of the FS can be naturally understood. 
Other satisfactory features include (a) the ratio of the 
spectral weight on the back side of the pocket to that on 
the main side, (b) the dispersion kink in the nodal direc- 
tion around « 0.05 eV, and (c) the shadow band disperses 
out as the energy is lowered below the Fermi energy as 
Fig. 4(b) shows because the shadow feature is "reflec- 
tion" of the main band with respect to the (0, tt) — (tt, 0) 
line. In the laser ARPES experiments by Meng et al. the 
shadow band was observed to disperse out as the binding 
energy increases. See the lower row of the plots b-d in 
the Fig. 1 of Ref£. 

Despite these satisfactory features of the dynamic fluc- 
tuations there are some discrepancies compared with ex- 
perimental observations. First of all, the current scenario 
requires quite long correlation length of order of £/a ~ 10 
for the Fermi arcs or Fermi pockets to appear. But, one 
of the present authors recently inverted the high resolu- 



tion laser ARPES from Bi2Sr2CaCu20s+5 in pseudogap 
state to extract the bosonic fluctuations spectrum shown 
in Fig. 6. It was found that the correlation length is of 
the order of £/a < O.l^ Although the Eliashberg func- 
tion was extracted in Bi2Sr2CaCu20s+(5 and the Fermi 
pocket/arc was observed in Bi2Sr2- 2; La 2 ;Cu06 , both ex- 
periments were carried out in the pseudogap state and 
this contradiction needs to be reconciled. 

Secondly, if the fluctuation spectrum of Fig. 6 is peaked 
at (tt, tt) with £/a ~ 10, then the transport \ tr ~ 2A * 3 
in the nodal direction because of the (1 — cosO) factor 
from the vertex correction. To our knowledge, this large 
Xtr was not observed in the resistivity measurements. 
Bok et al. concluded that the correlation length must 
be small, £ <C a, and the spectrum can not be from the 
(tt,tt). The enhancement of Xt r over A is not expected. 
An interesting point in this context though is the obser- 
vation by Schachinger and Carbotte^ They compared 
the a 2 F from infrared (IR) spectroscopy and ARPES 
and found that they agree well overall (after scaling) ex- 
cept that a 2 F from IR is larger than that from ARPES 
around 0.06 eV by the factor of approximately 2. This 
may be understood if the peak around 0.06 eV is dom- 
inantly from (tt, tt) and the rest of the spectrum is mo- 
mentum independent. However, this scenario seems to 
be at odds with the conclusion of Bok et al. 

As the coupling constant a increases, the electron 
Fermi surfaces disappear first leaving the hole Fermi sur- 
faces only as the two-pole approximation illustrates in 
Fig. 3. This topological change of the Fermi surface, 
however, was not obtained in the self-consistent calcula- 
tions because the iteration procedures failed to converge 
for a larger than approximately 0.22 eV. The a > 0.2 
eV and/or momentum dependent a is expected to give 
interesting results about the Fermi surface evolution as a 
function of doping, coupling constant, and temperature. 

It should be also interesting to check if one can under- 
stand the quantum oscillations under the applied mag- 
netic field with the current scenario. It is conceivable that 
the dynamically induced hole Fermi arcs/pockets are sup- 
pressed and the electron pockets are formed as the field 
is applied as the quantum oscillation experiments imply. 

Finally, we wish to note that Chang et al. also observed 
the back-side of the Fermi pocket in the pseudogap state 
in Lai.4gNdo.4Sro.i2Cu04 where the orthorhombic distor- 
tion is not the primary caused Recall that the previous 
observations of the shadow bands were found to be due 
to the orthorhombic structural distortion. 31 This struc- 
tural feature was separated out in Meng et al. Also the 
improved resolution of the laser ARPES facilitated their 
observation of the Fermi arcs and pockets. An interesting 
point is that the observed shadow band by Chang et al. 
was much stronger than Meng et al. and was more sym- 
metric with respect to the (tt, 0) — (0, tt) line. It remains to 
be sorted out what causes the differences between Meng 
et al. and Chang et al. 
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